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Abstract 

We study the radiative decay — * 7r°r/7 within the framework of a phe- 
nomenological approach in which the contributions of p-meson, chiral loop 
and ao-meson are considered. We analyze the interference effects between 
different contributions and utilizing the experimental branching ratio and in- 
variant ir°n mass spectrum for <f> — ► 7r r^7 decay we estimate the branching 
ratio of (ft — ► ao7 decay. 
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I. INTRODUCTION 



The nature and quark substructure of low mass scalar mesons, in particular those of 
isoscalar / (980) and isovector a (980), have been a subject of controversy in hadron spec- 
troscopy over the years. A variety of interpretations have been proposed for their structure, 
but, whether they are conventional qq states in quark model KK molecules [[|, or exotic 
multiquark q 2 q 2 states || have not been established yet. 

It was suggested by Achasov and Ivanchenko 0] that the radiative decays of meson to 
pseudoscalar mesons, — > 7r°7r°7 and — > 7r°?77, offer the possibility of obtaining informa- 
tion on the nature of /o(980) and ao(980) mesons, respectively. Close, Isgur and Kumano f| 
noted that the radiative decays — > £7, where S = fo or clq, can be utilized to differentiate 
among various models of the structure of these scalar mesons. They shown that although 
the transition rates r(0 — > /07) and T(<f> — > 007) depends on the unknown dynamics, the 
ratio of the decay rates T((f> — > /o7)/r(0 — > a 7) will be sensitive to the spatial distribution 
of quarks and the spatial wavefunctions of scalar mesons, and thus it provides an experi- 
mental test which distinguishes between alternative explanations of their structure. In both 
the KK molecule, and four-quark cluster q 2 q 2 pictures of the structure of scalar mesons, 
it is now generally accepted that the radiative decays — - > 5*7 proceed by the mechanism 
in which and S both couple to an intermediate K + K~ loop in the chain of reactions 
— > K + K~ — > K + K^ r y — > S7 The corresponding decay rates were obtained as || 

BR((f)^ / (980)7) = BR(<p -> a (980)7) ~ (2.0 ±0.5) x x F 2 (R) , (1) 

where F 2 (R) is a form factor that depends on the spatial wavefunctions of the scalar mesons, 
and F 2 (R) = 1 in point-like effective field theory calculations. In obtaining this result it is 
further assumed that both / (980) and a (980) are isospin eigenstates, and gg K+K -/4:n = 
0.58 GeV 2 . If /o(980) and ao(980) are spatially extended KK molecules with R > Aqq D , 
then F 2 (R) < 1 and the resulting branching ratios are BR(0 — > S7) ~ (0.4 — 1) x 10 -4 0. 
If the meson states are compact four-quark clusters confined within R ~ ^qcd in q 2 q 2 
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picture, then the resulting branching ratios are expected as BR(0 — > Sj) ~ 2 x 10~ 4 ||||. In 
addition to the predictions about the absolute branching ratios of — > S7 decays, there is 
also the weaker result that follows from Eq. 1, that is theoretically BR(<f) — > /07) /BR((f) — > 
a 7) = 1. 

The very accurate data that have been obtained from Novosibirsk SND |],[7] and CMD-2 
H Collaborations give the following branching ratios for — > 7r°7r°7 and — > ir r)^f decays: 
BR(0 -> 7r°7r°7) = (1.221 ± 0.098 ± 0.061) x 10~ 4 §, BR(0 -> tt ^) = (0.88 ± 0.14 ± 
0.09) x 10~ 4 0, and BR(0 -> vr 7r°7) = (0.92 ± 0.08 ± 0.06) x 10~ 4 , BR(0 -> ttV) = 
(0.90 ± 0.24 ± 0.10) x 10~ 4 |J, where the first error is statistical and the second one is 
systematic. From the analysis of their experimental results, Achasov et al. |J concluded 
that the / 7 mechanism dominates the — > 7r°7r°7 decay and the contributions coming from 
o"7 and p°7r° intermediate states are small. Neglecting these contributions they obtained 
the branching ratio BR(0 — > /07) = (3.5 ± 0.3io;i) x 10~ 4 , where the second errors are 
systematic and includes uncertainty in the interference terms. This branching ratio was also 
obtained by Akhmetshin et al. from their data of <fi — > 7r°7r°7 decay as BR(0 — > f (980)7) = 
(2.90 ±0.21 ± 1.54) x 10~ 4 §]. The <p — > / (980)7 decay branching ratio is now quoted as an 
average of these two results which is BR(0 — > / (980)7) = (3.4 ± 0.4) x 10 -4 0. Achasov 
et al., on the other hand, by assuming that ao7 intermediate state dominates the — > 7r°?77 
decay, and the contributions from other decay mechanisms, for example <fi — > p°7r°, p° — > 177, 
can be neglected, obtained the value for the branching ratio of the decay — > a (980)7 as 
BR(0 — > ao(980)7) = (0.88 ±0.17) x 10~ 4 from the data of their — > 7r ?77 decay experiment 
|8j. Therefore, these experimental results suggest for the ratio of the decay rates of the decays 
— > /o(980)7 and — >■ a (980)7 the experimental value T(0 — > f j)/T((f> — > 007) = 3.8±1.2, 
with which the theoretical result T(0 — > / O 7)/r(0 — > a 7) — 1 is not in agreement. 

In this paper, we study the role of the scalar ao(980) meson in the mechanism of the 
radiative — > 7r ^7 decay employing a phenomenological framework, in which we try to 
assess the roles of different processes and the contributions to the decay rate coming from 
their amplitude in the mechanism of this decay, and this way we attempt to estimate the 
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branching ratio BR(0 — > a® (980)7) °f the decay — > ao (980)7 horn the experimental data of 
the radiative <ft tt°^7 decay. In our analysis, we use the experimental data of Novosibirsk 
SND Collaboration 0. 

II. FORMALISM 

The radiative decay process of the type V° — > P°P°7 where V and P belong to the 
lowest multiplets of vector (V) and pseudoscalar (P) mesons have been studied using differ- 
ent approaches. In particular, the branching ratio BR(<p — > 7i i] , ~f) has been calculated as 
BR((j) — > 7T T]j)vdm — 5.4 x 10" 6 |I(J by considering intermediate vector meson contribu- 
tion only, and as BR((p — > 7r ^7) x = 3 x 10~ 5 [0] by using chiral loop model within the 
framework of chiral phenomenological Lagrangians. 

In order to include the scalar ao meson resonance pole in the decay mechanism for the 
calculation of the decay rate of <p ~~ > 7r°?77 decay in a phenomenological framework we 
consider two different approaches. In both approaches, we do not make any assumptions 
about the structure of the scalar ao meson. 

In the first approach, which we name Model I, we include the scalar ao resonance in an 
ad hoc manner. In this approach, we assume that the mechanism of the — > ir -)]^ decay 
consists of the reactions shown by the diagrams in Fig. 1. We describe the </>a 7 and a ir°r) 
vertices in Fig. 1(c) by the phenomenological Lagrangians 

£*o<>7 = TF^o7[^ Q< / ,/3 ^^/3 - d a (f) l3 df 3 A a \a (2) 

and 

respectively, which also serve to define the coupling constants a^, ao7 and g ao7rrr Since there 
are no direct experimental results relating to 0a O 7-vertex, but only an upper limit for the 
branching ratio of the decay — * a 7 is quoted in Review of Particle Properties as BR(<fi — > 
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ao7) < 5 x 1CT 3 0, we determine this coupling constant in our calculation by employing 
the experimental branching ratio of the <fi — ► 7r°?77 decay. The decay rates for the <fi — > a 7 
and the ao — > tt°i] decays resulting from the above Lagrangians are 



T(<j) -> a 7) 



24vr 



MI 



PQ07 ' 



(4) 



and 



T(a -> vr°?7) 



M 2 

a 



l _ (M^o - 



M 2 

a 



(5) 



167rM ao ^ [ 

respectively. We use the value r ao = (0.069 ±0.011) GeV determined by E852 collaboration 
at BNL Hl2] , and we obtain the coupling constant g ao -nr) as g ao -wr) = (2.32 ± 0.18) GeV. 

The (ppir vertex in Fig. 1(a) is conventionally described by the phenomenological La- 
grangian 



(6) 



and the coupling constant <^ p7r was determined by Achasov and Gubin using the data on 
the decay (f) — > pit — > 7r + 7r~7r° || as g<p p7T = (0.811 ± 0.081) GeV~ x |13]. The pir^ vertex in 
Fig. 1(a) is described by the phenomenological Lagrangian 



(7) 



and the coupling constant g pm is then obtained from the experimental partial width of the 
radiative decay p — > 777 || as g pr]1 = (0.45 ± 0.07) GeV^ 1 . 

We describe the (j)KK vertex in Fig. 1(b) by the phenomenological Lagrangian 



-<t>K+K- 



(8) 



which results from the standard chiral Lagrangians in the lowest order of chiral perturbation 
theory flH |. The decay rate for the <fi — > K + K~ decay resulting from this Lagrangian is 

3/2 

M, 1 - I 1 i • (9) 



K + K~ 
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We utilize the experimental value for the branching ratio BR(<f) — > K + K~) = (0.492±0.007) 
for the decay 4> — > K + K~ ||, and determine the coupling constant g^KK as g^KK = (4.59 ± 
0.05). For the four pseudoscalar KKtit] amplitude, we use the result obtained in standard 
chiral perturbation theory |15| which is 

M{K + K~ -> 7r°T]) = (m*o v - ^ M K^j • (10) 

where rj — r{ mixing is neglected, M n o v is the invariant mass of the 7r°r] system, and we use 
f n = 92.4 MeV. We therefore obtain the amplitude for the diagram in Fig. 1(b) as 



M 



e g<t>KK 



[ip ■ k) (e • u) - (p ■ e) (k ■ u)] I (a, b)M (K + K~ -> n r]) 



(11) 



where (u,p) and (e, k) are the polarizations and four-momenta of the meson and the 
photon, respectively, and a = M^/M^, b = M 2 0ri /M^ with the invariant mass of the final 
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] r] system given by M^o = (qi + q-i) 2 = (p — k) 2 . The loop function is defined as fl6| , |i~7 
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On the other hand, the introduction of the ao amplitude as in Fig. 1(c) may be considered 
not to be very realistic since this diagram implies direct quark transition and it thus makes 
a very small contribution because of OZI supression. Indeed, it has been shown that the 
scalar resonances /o(980) and ao(980) can be excited from the chiral loops, with the loop 
iteration provided by the Bethe-Salpeter equation using a kernel from the lowest order chiral 



Lagrangian p8| . Furthermore, it has been argued that the experimental data obtained in 
Novosibirsk give resonable arguments in favour of the one- loop mechanism for — ► K + K~ — > 
ao7 and — > K + K~ — > / 7 decays [0. Therefore, in our second approach to study the 
radiative — > ir ^ decay, which we name Model II, we assume that this decay proceeds 
through the reactions — > p°7i° — > 7r°?77, — > K + K~^ — > 7r°?77, and — > ao7 — > 7r r)~f 
where the last reaction proceeds by a two-step mechanisms with ao coupling to with 
intermediate KK states. We show the processes contributing to the — > 7r°?77 decay 
amplitude diagramatically in Fig. 2. We note that we do not make any assumptions about 
the structure of the ao meson, and only assume that the and ao mesons both couple 
strongly to the K + K~ system, as a result of which there is an amplitude for the decay 
- ► o 7 to proceed through the charged kaon loop independent of the nature and the 
dynamical structure of a meson. The only new ingredient required in our second approach 
is the K + K~ ao vertex which we assume is described by the phenomenological Lagrangian 

C ao K+K- = 9a K+K~ K + K~a . (14) 
The decay width of ao that follows from this Lagrangian is 



T(a -> K + K~) - 9 ^ K+K ~ 



'2M K \ 2 



M ao 



1/2 

(15) 
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which is usually considered to define the coupling constant g au K+K- ■ In our second approach, 
we calculate the decay rate of — > 7r°?77 decay using the diagrams shown Fig. 2, and by 
utilizing the experimental value of this decay rate determine the coupling constant g ao K+K- ■ 
In our calculation of the invariant amplitudes we make the replacement p 2 — M 2 — > p 2 — 
M 2 + iMT in ao and p° propagators and use the experimental value T p = (150.2 ±0.8) MeV 
for p° meson, because using a g 2 -dependent width did not affect our results appreciably. 
However, since the mass of the M ao meson is very close to the K + K~ threshold, this induces 
a strong energy dependence in the width of the ao-meson. We follow the widely accepted 
option to deal with this problem that was proposed by Flatte PDA based on a coupled channel 



(nr), KK) description of the ao resonance, and parametrize the ao width as 



= O^g 2 ) 0{yff- (M„ +M V )) 



+ig K ^Ml - 6(2M K - yj<?) + g K ^q 2 /4 - M\ 6{^q* - 2M K ) , (16) 

where 



r "> 2 ) = 1fi !w /a v / fe 2 -( M ^+ M .) 2 ]fe 2 -( M ^- M .) 2 ] • ( 17 ) 
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We use the Flatte parameter g K g as = (0.22 ± 0.04) which was determined by E852 
Collaboration at BNL by a fit to data in their experiment in which they determined the 
parameters of the do meson [|T^]. In all our calculations and in the analysis of the experi- 
mental invariant ir r] mass spectrum of the <fi — > 7r°?77 decay, therefore, for the mass of ao 
meson we use the value M ao = (0.991 ± 0.0025) GeV as determined in this experiment by 
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We express the invariant amplitude .M(E 7 , Ex) for the decay <fi — ► 7r°?77 as M. = M. a + 
Mb + Mc in Model I and as M = M' a + M' b + M' c in Model II, where M a , Mb, and M c are 
the invariant amplitudes resulting from the diagrams (a), (b), and (c) in Fig. 1 respectively, 
describing the Model I, and M' a , M' b and M' c are the invariant amplitudes corresponding 
to the diagrams (a), (b) and (c) in Fig. 2, respectively, defining the Model II. This way we 
take the interference between different reactions contributing to the decay <fi — > 7r°?77 into 
account. Then the differential decay probability for <p — > 7r°?77 for an unpolarized 4> meson 
at rest is given as 

dT = _i L I M I 2 (18) 

dE 1 dE l (2tt) 3 8M 1 1 ' K } 

where E 7 and E x are the photon and pion energies respectively. We perform an average over 
the spin states of meson and a sum over the polarization states of the photon. The decay 
width r(0 — > TT i]'~f) is then obtained by integration 

T = / dE 1 / dE 1 —— r (19) 

where the minimum photon energy is E 7)m j n . = and the maximum photon energy is given 
as E ljmax . = [M? — (M n o + Mj 1 ) 2 ]/2M l f ) . The maximum and minimum values for pion energy 
Ex are given by 



—-L-—[-2E%M4 - M^Ml + Ml, - Ml) + £ 7 (3Mj + M 2 n0 - M*) 
±E^4E*M* + Ml + (Mj - M*)* - 2Mj(Mj + Mj) + \E-M (~M;, + M^ + M[) 



III. RESULTS AND DISCUSSION 

In order to determine the coupling constants g^-y m Model I and g ao K+K- in Model II, 
we use the experimental value of the branching ratio for the radiative decay <fi — > 7r°?77 H in 
our calculation of this decay rate, and this way we obtain a quadric equation for the coupling 
constant g^ ao -y in Model I and another quadric equation for the coupling constant g ao K+K- 
in Model II. In these quadric equations the coefficient of quadric term results from the ao 
meson amplitude contribution, and the coefficient of the linear term from the interference 
of the ao meson amplitude with the vector meson dominance and the loop amplitudes. We 
then predict and study the invariant mass distribution dBR/dM n o v = (M 7T o v /M ( f ) )dBR/dE.y 
for the radiative decay — > 7r°?77 in our phenomenological approach using the values of 
coupling constants g^o-y in Model I and g ao K+K- in Model II that we obtain and compare 
our results with the experimental invariant 7r°7] mass spectrum for the decay <fi — > 7r r]7 [[?[]. 

In Model I, we obtain for the coupling constant 5^07 the values g^ao-y = (0.24 ± 0.06) 
and g^ao-y = ( — 1-3 ± 0.3). In Fig. 3 we plot the distribution dBR/dM 7T o v choosing the 
coupling constant g0 ao7 = (0.24 ± 0.06) in which we also indicate the contributions coming 
from the different reactions shown diagramatically in Fig. 1, as well as the contribution 
of the total amplitude which includes the interference term as well. Our Model gives a 
reasonable prediction for the spectrum over most of the range of the invariant mass M n o v 
except in its higher part where the expected enhancement due to the contribution of ao 
resonance is not produced. The distribution dBR/dM n o v we obtain for the other root, 
that is for = (—1.3 ± 0.3), is even poorer which we do not show. We obtain the 

branching ratio for the decay — > ao7 using for the coupling constant 0^07 — (0.24 ± 0.06) 
as BR((p — > 007) = (0.2 ± 0.1) x 1CT 5 . However, because of the fact that Model I does 
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not produce a satisfactory description of the experimental invariant M n o„ mass spectrum 



for the decay ~> 7i"°^7, we cannot consider the value of the coupling constant 



ao7 



(0.24 ± 0.06) and the resulting branching ratio BR(<p — > a 7) = (0.2 ± 0.1) x 10~ 5 too 
seriously. In this connection, we like to note that the coupling constant g</, ao7 has been 
calculated [^] employing QCD sum rules and utilizing c<j0-mixing by studying the three 
point 0a O 7-correlation function as 5^07 = (0.11 ±0.03). Furthermore, the photoproduction 
cross section of p° mesons on photon targets near threshold has been shown |22| to be mainly 
given by cr-exchange, and assuming vector meson dominance of the electromagnetic current 
the value of the coupling constant g pal = 2.71 was deduced. In the study of the structure of 
the meson photoproduction amplitude on nucleons near threshold based on the one-meson 
exchange and Pomeron-exchange mechanism ||23|| , this value of the coupling constant g p<T1 



was used to calculate the coupling constant 5^07 by invoking unitary symmetry arguments 
as I 9<t>a Q i |= 0.16 by assuming that a, fo, and ao are numbers of a unitary nonet, which 
is not without problems. In the light of our discussion of Model I and the resulting poor 
invariant mass spectrum of the decay — > tt -!]^, it may be argued that all these values for 
the coupling constant g^ao-y should be considered with some caution. 

We follow the same procedure in Model II, and obtain the coupling constant g ao K+K- 
as g ao K+K- = (—1-5 ± 0.3) GeV and g ao K+K- = (3.0 ± 0.4) GeV utilizing the experimental 
value of the ~> decay rate . We then plot the resulting invariant mass distribu- 

tion for the decay — > tt ^ and compare it with the experimental result [0]. In Fig. 
4 and in Fig. 5 we plot the distribution dBR/dM n o v for the radiative decay — > 7r r)~f 
in our phenomenological approach choosing coupling constant g ao K+K- = —1-5 GeV and 
9o, k+k- — 3.0 GeV, respectively, as a function of the invariant mass M^o^ of the 7r°r) sys- 
tem. In these figures, as before, we also indicate the contributions coming from the different 
reactions <fi —>■ p°7r° — > 7r°?77, — > K + K^ r y — > 7r°?77, and — > a 7 — > TT°i]'-f in Model II 
shown diagramatically in Fig. 2, as well as the contribution of the total amplitude which 
includes the interference terms as well. We note that for g ao K+K- = —1-5 GeV our predic- 
tion for the invariant mass spectrum is in good agreement with the experimental result and 



not only the overall shape and feature but also the enhancement due to the contribution of 
the do resonance is well produced. Therefore, from the analysis of the spectrum obtained 
with the coupling constants g ao K+K- — —1-5 GeV and g ao K+K- — 3.0 GeV in Fig. 4 and 
5, respectively, we may decide in favour of the value g m K+K- = —1-5 GeV, and we may 
state that the experimental data within the framework of our phenomenological approach 
in Model II, suggest for the coupling constant g ao K+K- the value g ao K+K- = (—1-5 ± 0.3) 
GeV. Furthermore, we note that Model II provides a better way as compared to Model I in 
order to incorporate the do meson into the mechanism of the <fi — > 7r r]7 decay, and thus our 
result gives further support for the approach in which a® meson state arises as a dynamical 
state and for the two-step mechanism for its coupling to <fi meson with intermediate KK 
states. In a previous work |24|, we studied the radiative <fi — > 7r°7r°7 decay also within the 



framework of a phenomenological approach in which we considered the contributions of a 
meson, p meson and /o(980) meson. We analyzed the interference effects between different 
contributions using the experimental results of SND Collaboration. In that analysis, the 
coupling constant gj K+K- appeared with positive sign, that is gf K+K- > 0. In the present 
analysis, we obtained the coupling constant g ao K+K- as g ao K+K- < 0- We like to note that 
this result about the relative phase between gf K+K- and g ao K+K- is consistent with the 
result obtained in the q 2 q 2 model where gf x+K- = —g ao K+K- H,i,Pl • 

We then calculate the decay rate T(<f) — > a 7) by assuming the two-step one-loop mech- 
anism — > K + K~ — > ao7 for the decay <fi — *> ao7- We show this mechanism diagramatically 
in Fig. 6. The decay rate that follows from these diagrams is given by 

r(» -. a„7) = a9l %^:" KK ^ I (» - Wa,b) I 2 . (20) 

where a = Mj/Mj,, b = M 2 o /Mj,, and the loop function I(a, b) is defined in Eq. 11. We 
use the coupling constant g ao KK = (—1-5 ± 0.3) GeV and then obtain the decay rate for 
T(0 — > aoj) decay as T(<j) — > ao7) = (0.51 ± 0.09) KeV and the branching ratio as BR{<p — > 
a ol) — (1-1 ±0.2) x 10~ 4 . This result when combined with the experimental value BR(<f) 
/ o7 ) = (3.4 ±0.4) x 10~ 4 gives the ratio T(0 f Q i)/Y{<}> -> a 7) = (3.0 ±0.9), which is still 
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larger than the theoretical value of l.The value BR(<p — > ao7) = (1.1 ±0.2) x 10~ 4 we obtain 
for the branching ratio of the decay ~> a oli although not in disagreement, is somewhat 
larger than the result of Achasov et al. which is BR(<fi — > a^) = (0.88 ± 0.17) x 10~ 4 [0. 
Consequently, we obtain a smaller value for the experimental ratio r(0 — > /o7)/T(0 — > ao7) 
than the previous result which is r(0 — > f j)/T((j) — > a 7) = (3.8 ± 1.2). We like to note, 
however, that Achasov et al [0] assumed that ao7 intermediate state dominates the — > 7r°r)~f 
decay and they neglected the contributions coming from other intermediate states in their 
analysis. On the other hand, in our analysis we include the contributions coming not only 
from ao7 intermediate state but also from p°7r° and K + K~ intermediate states as well as 
from their interference. 

In order to obtain the branching ratio for — > a 7 decay utilizing the experimental 
value of the decay rate and the experimental invariant ir°r] mass spectrum for the — > 7r r)~f 
decay, we use a phenomenological approach. In our analysis, we employ point-like effective 
field theory and we deduce from the experimental data the coupling constant of point-like 
a Q meson coupled to point-like K meson. The coupling constant that we obtain g ao K+K- — 
(-1.5 ± 0.3) GeV results in the branching ratio BR(<p -> a 7) = (1.1 ± 0.2) x 10~ 4 for 
— > a 7 decay using one loop — > K + K~ a 7 mechanism. We can, therefore, suggest 
that the branching ratio for the decay — > ao7 that is used in the literature should be 
revised as BR(<p —>■ a^) = (1.1 ± 0.2) x 10~ 4 x F 2 (R). We can also assert that our analysis 
suggests a lower value (3.0 ±0.9) for the ratio T(0 — > /o7)/r(0 — > a 7) than used previously. 

In a very recent paper |26j the KLOE Collaboration presented the result of their exper- 
imental study of — > 7r°^7 decay. In their analysis, they conclude that — > tt ^ decay 
is dominated by the process — > ao7 and from a fit to 7T°r) invariant mass spectrum they 
find the branching ratio BR(<fi — > a 7) = (7.4 ± 0.7) x 10 -5 for the — > a 7 decay. In 
view of this finding, our result BR(4> — > a ^) = (1.1 ± 0.2) x 10 -4 x F 2 (R) implies that 
F 2 (R) = (0.7 ± 0.2). This value supports the view that the structure of a meson is a 
combination of a KK molecule with a compact q 2 q 2 core 
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FIGURES 
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FIG. 1. Diagrams for the decay (f> — > tt 7]j in Model I. 
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FIG. 2. Diagrams for the decay <j) — > 7r 777 in Model II. 
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FIG. 3. The n rj invariant mass spectrum for the decay <p —> tt 777 in Model I. The contributions 
of different terms are indicated. 
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FIG. 4. The ir r/ invariant mass spectrum for the decay 4> —> n nl f° r ZoqK+k- = — 1-5 GeV in 
Model II. The contributions of different terms are indicated. 
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FIG. 5. The it rj invariant mass spectrum for the decay (j) — ► 7r°777 for g ao i<-+^- = 3.0 GeV in 
Model II. The contributions of different terms are indicated. 
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FIG. 6. Diagrams for the decay (p —> ^07- 
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